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Abstract 

We review the recent experimental and theoretical progress in the determination of \Vud\ 
and |Eu<i|, and the status of the most stringent test of CKM unitarity. Future prospects 
on \Vcd\ and \Vcs\ are also briefly discussed. 
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1 Introduction 


The most precise constraints on the size of the elements of the CKM matrix [1, 2] are 
extracted from the low-energy s —>■ m and d —>■ m semileptonic transitions. Combining the 
precise determinations of \Vud\ and | 14 s| extracted from these processes we can perform 
the most stringent test of CKM unitarity, namely we can probe the validity of the relation 

\Vud? + \Vus? + \Vub\^ = l. ( 1 ) 

at the 0.1% level. In particular, the best determination of \Vus\ is obtained from semilep¬ 
tonic K decays and K 12 ), while the most stringent constraints on \ Vud\ are obtained 
from superallowed Fermi transitions (and, to a minor extent, from neutron and pion beta 
decays). As we will discuss in the following, in the last few years there has been a substan¬ 
tial progress - both from the theoretical and the experimental side - in the determination 
of these two matrix elements. 

In all cases the key observation which allow a precise extraction of the CKM factors 
is the conservation of the vector current at zero momentum transfer in the SU{N) limit 
and the non-renormalization theorem. The latter implies that the relevant hadronic form 
factors are completely determined up to tiny second order isospin-breaking corrections in 
the d —>• M case [3] or SU(3)-breaking corrections in the s —> m case [4]. As a result of 
this fortunate situation, the accuracy on | 14 s| is approaching the 1% level and the one on 
\Vud\ is already below the 0.05% level. 

The present accuracy on \Vud\ and |14s| is such that the contribution of |14fep ~ 2x10“® 
in the relation m can safely be neglected, and the uncertainty of the first two terms is 
comparable. In other words, to a high degree of accuracy we can set 

Vus sin 0^ , ^ud cos 0^ , (2) 

as in the original Cabibbo theory [1], and \Vus\ and \Vud\ provide two independent deter¬ 
minations of the Cabibbo angle both around the 1% level. 

In the following sections we will review the determinations of |14s| and \Vud\ from the 
main observables mentioned above. We will also briefly analyze alternative strategies to 
extract |14s| from r and hyperon decays, as well as the future prospects on iWsI and 
\Vcd\- The main recent results on |14s| and \Vud\ are then summarized and combined in 
the last section, where we will discuss the accuracy to which Eq. m is satisfied and we 
will provide a final global estimate of the Cabibbo angle. 

2 The extraction of Vud 

The value of Vud has been extracted from 1) super-allowed, 0'*' —> O’*", nuclear beta decays, 
2) neutron beta decays, n —>• peu, and 3) pion beta decay 7r+ —>• The latter two, 

subsequently discussed in this report, have smaller overall theoretical uncertainties and 
may in the long term be better ways to obtain Vud] hut currently, only super-allowed beta 
decays determine Vud to better than 0.05%. 


2 


2.1 Super-allowed Fermi transitions 

The so-called super-allowed, O’*" —*■ 0^, Fermi transitions between nuclei are very spe¬ 
cial [5]. Because they proceed (at the tree level) through pure weak vector current inter¬ 
actions, which are conserved in the limit; they are not renormalized by strong 

interactions at = 0. Hence, they are ideally suited for cleanly extracting V^d with high 
precision. Corrections due to 7 ^ 0 and rud 7 ^ are negligibly small; so, one needs 
only to control uncertainties in the electroweak radiative corrections, isospin violating 
electromagnetic effects and nuclear structure dependence. How well that can be done is 
the subject of this section. 

Last year, the prevailing value of Vud obtained by averaging the nine best measured 
super-allowed /3-decays was [6,7] 

Kd = 0.9740(1)(3)(4) ^ 0.9740(5) (2004 value) (3) 


where the errors are experimental, nuclear theory and quantum loop corrections. The very 
small experimental error illustrates the power of this averaging procedure. The largest 
uncertainty, associated with weak axial-vector induced loop effects [ 8 ] primarily through 
7 IT box diagrams, represents model dependent hadronic effects which until recently [9] 
were thought to be essentially irreducible or at least very difficult to reduce. 

Two developments have led to a recent improvement in Vud by nearly a factor of 2. 
First, new global studies of super-allowed /3-decays by Hardy and Towner [10], and by 
Savard et al. [11], have provided a more consistent treatment of Q values and lifetimes 
used in ft determinations, which in turn give Vud via the master formula 

2984.48(5) s 
ft{l + RC) 

In that expression, RC designates the total effect of all radiative corrections from quantum 
loops as well as nuclear structure and isospin violating effects. RC is nucleus dependent, 
ranging from about -|-3.1% to -|-3.6% for the nine best measured super-allowed decays. 
That difference is of critical importance in bringing the values of Vud obtained from sepa¬ 
rate decays into agreement with one another. The magnitude of the corrections is essential 
for establishing unitarity, as we shall see. 

A second major advance in the determination of Vud stems from a new study of the 
quantum loop corrections coming from the problematic 7 W box diagram due to weak 
axial-vector contributions. Previously, those effects, along with other smaller axial-vector 
current contributions, were found to shift the RC by about 

( 5 ) 

where Ag = —0.34 is a one-loop QCD correction to the short-distance logarithmic loop 
contribution and Csorn — 0.8gA(/^p + /^n) — 0.9 represents long-distance loop effects. The 
problematic intermediate loop momentum region was roughly estimated by employing 


a 


In-h Ag -\- 2C'Born 

rriA 



( 4 ) 
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rriA — 1.2GeV in the log, while the crudely obtained error of ±8 x 10““^ in that quantity 
(which leads to ±4 x 10“^ in Vud) was found [6,8] by allowing the cut-off scale to vary 
up or down by a factor of 2. 

A new analysis [9] of the yW box diagram now divides the loop momentum into 3 
integration regions: 


(1.5 GeV)2 <Qj<oo 
(0.82 GeV)2 < Qjj < (1.5 GeV)^ 

0 < Qhi < (0.82 GeVf 

The evaluation of region I has been supplemented by 3-loop QGD corrections to the 
leading term in the short-distance operator product expansion, rendering it effectively 
error free and, more important, allowing a smooth extrapolation to lower Q^. Region 
II has been evaluated using interpolating vector and axial-vector resonances, a proce¬ 
dure motivated by large Nc QGD and vector meson dominance. That prescription has 
been well tested in other calculations; nevertheless, a conservative ±100% uncertainty 
has been assigned to that part of the calculation. Finally, region III was evaluated using 
well-measured nucleon dipole form factors and assigned a ±10% uncertainty. Those im¬ 
provements have reduced the theoretical quantum loop uncertainty in Vud from a crude 
±4 X 10“'^ to a more defensible ±1.9 x 10“'^, about a factor of 2 improvement. Further 
error reduction may be possible if future lattice calculations can confirm the interpolating 
resonance approach, since the uncertainty from intermediate momenta is still dominant. 

The overall shift in Vud due to the new evaluation of radiative corrections is relatively 
small. However, the error reduction is more signihcant. Updating the most recent ft 
values [11] with the new RG results leads to the Vud values given in Table 1. Gombining 
all errors in quadrature now gives the weighted average [9] 

Vud = 0.97377(27) (2005 value) (6) 

The central value has not shifted very much [see Eq. (7)], but the error has been reduced 
by nearly a factor of 2. 

So far the situation for Vud looks very good; however, we caution that the recent re¬ 
measurement of the Q value for [11] has implied a substantial (~ 2a) shift in the 
corresponding extraction of Vud [from 0.97363(50) to 0.97280(43)]. As a result, the overall 
consistency of the Vud values extracted from the various nuclei is not as good as it was 
about one year ago. This fact could be interpreted as an indication of possible problems 
with the Z-dependent radiative corrections or the Q values of other superallowed decays. 
The shift in the total average of \Vud\ implied by the new data is not significant, but 
changes in the other superallowed Q values could have more substantial effect. 

The superallowed beta decays have now reached the very impressive ±0.03% level 
of precision in their determination of Vud- Further studies of those reactions are clearly 
warranted, both to reduce the error and to clarify the new anomaly. In addition, 
future high statistics neutron-studies [6,12] of and qa may be able to reach a level of 
precision for Vud comparable to Eq. but without the nuclear physics complications. 
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Nucleus 

ft (sec) 

Vud 

IOC 

3039.5(47) 

0.97381(77)(15)(19) 

140 

3043.3(19) 

0.97368(39)(15)(19) 

26 Al 

3036.8(11) 

0.97406(23)(15)(19) 

34C1 

3050.0(12) 

0.97412(26)(15)(19) 

38k 

3051.1(10) 

0.97404(26)(15)(19) 

42 Sc 

3046.8(12) 

0.97330(32)(15)(19) 

46v 

3050.7(12) 

0.97280(34)(15)(19) 

soMn 

3045.8(16) 

0.97367(41)(15)(19) 

54Co 

3048.4(11) 

0.97373(40)(15)(19) 

weighted ave. 

0.97377(11)(15)(19) 


Table 1: Values of Vud implied by various precisely measured superallowed nuclear beta 
decays [9]. The ft values are taken from a recent update by Savard et al. [11], Uncertain¬ 
ties in Vud correspond to 1) nuclear structure and Z‘^o? uncertainties added in quadrature 
with the ft error, 2) a common error assigned to nuclear coulomb distortion effects, and 
3) a recently reduced (common) uncertainty in the radiative corrections from quantum 
loop effects. Only the hrst error is used to obtain the weighted average. 

Those measurements, which are difficult but well worth the effort, will be discussed in the 
next subsection. 

2.2 Neutron /3-decay 

The beta decay of the neutron allows a determination of \Vud\ free from the nuclear 
structure effects of superallowed beta decays. Fixing the Fermi constant from the muon 
decay, within the Standard Model (SM) we can describe the neutron /9-decay in terms of 
two parameters. One of them is \Vud\-i the other is the ratio of axial and vector coupling 
constants relevant to the n — > pev transition: A = Qa/Q v- The determination of \Vud\ is 
then based on two experimental inputs: the neutron lifetime and A. 

The neutron lifetime can be written as 


" |KdPGMl + 3A2)(l + A^)/^ ’ ' ' 

where = 1.71335(15) [6,13] is the phase space factor (corrected for the model inde¬ 
pendent radiative corrections), = 0.0239(4) denotes the model dependent radiative 
corrections to the neutron decay rate [9,14] and K is an appropriate normalization con¬ 
stant (see e.g. Ref. [15]). 

The most precise experimental information on A is derived from the /3-asymmetry 
coefficient Aq, which describes the correlation between the neutron spin and the electron 
momentum. To a minor extent, also the correlation between neutrino and electron mo¬ 
menta and the correlation between neutron spin and proton momentum are sensitive to 
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A. The neutron f3 decay is therefore an overconstrained system. In principle, A could 
also be determined from lattice QCD; however, at present the results of the most precise 
calculations are affected by 0(30%) errors. 

Since the overall uncertainty in \Vud\ is dominated by the experimental errors on 
and Ao, in the following we restrict our discussion on these two main observables. 


2.2.1 First input: A from the /3-asymmetry coefficient Aq 

The coefficient Aq is linked to the probability that an electron is emitted with angle ■d 
with respect to the neutron spin polarization P: 

lF(d) = 1 +-PAlocos(d) (8) 

c 

where n/c is the electron velocity expressed in fractions of the speed of light. Neglecting 
order 1% corrections, Aq is a simple function of A: 


Aq — 


o A(A + 1) 

1 + 3A2 ’ 


(9) 


where we have assumed that A is real. 

The most precise measurement of Aq, recently obtained by means of the instrument 
PERKEO, results in A = —1.2739(19) [17]. A recent repetition of this experiment confirms 
this value. In this experiment, the total correction to the raw data is 0.4% and the error 
contribution to \Vud\ is ±0.0007 [18]. Earlier experiments [19-21] gave significantly lower 
values for Aq. Averaging over recent experiments using polarizations of more than 90%, 
the Particle Data Group [22] obtains the value A = —1.2720(18). 

About half a dozen new instruments have been planned or are under construction 
for measurements of Aq and the other /3-asymmetry coefficients at the sub-10“^ level. 
Once this program will be completed, the error on \Vud\ due to the determination of A 
will be subleading with respect to the theoretical uncertainties of radiative corrections. 
Better neutron sources, in particular for high fluxes of cold and high densities of ultra-cold 
neutrons will boost the fundamental studies in this field. Major improvements both in 
neutron flux and degree of neutron polarization have already been made: first, a ballistic 
supermirror guide at the Institute-Lane Langevin in Grenoble gives an increase of about 
a factor of 4 in the cold neutron flux [23]. Second, a new arrangement of two supermirror 
polarizers allows to achieve an unprecedented degree of neutron polarization P of between 
99.5% and 100% over the full cross section of the beam [24]. Future trends have been 
presented at the workshop “Quark-mixing, GKM Unitarity” [12] and the “International 
Gonference on Precision Measurements with Slow Neutrons” [25]. 


2.2.2 Second input: the neutron lifetime 

The world average value for the neutron mean lifetime includes about a dozen individual 
measurements, using different techniques summarized as “in-beam methods” and “bottle 
methods”. The “in-beam methods” count the neutron decay product near a slow neutron 
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beam whereas the “bottle methods” store ultra-cold neutrons and count the neutrons that 


survived a certain time interval. One of the key points of neutron storage experiments is 
the control of neutron losses: in all bottle experiments performed so far, corrections of 
several percent are necessary because of the losses of ultra-cold neutrons at the material- 
trap walls. The best beam experiments have overall larger errors, but they need to apply 
smaller corrections. 

In 2004, all measurements agreed nicely with a x^/d.o.f.=0.95. The world average 
value was dominated by a single experiment [26] and the value adopted by the Particle 
Data Group was = 885.7(8)s [22]. The error contribution to \ Vud\ from this value of the 
neutron lifetime is ±0.0004, which is subleading with respect to the error induced by A. 

The situation has changed a few months ago, after the results of a new bottle experi¬ 
ment where the storage time for ultra-cold neutrons is very close to the neutron lifetime. 
The result reported by this experiment is = 878.5(8)s [27], which differs by more than 
6 cr from the Particle Data Group average. Such a change in the neutron lifetime would 
have a very signihcant effect on GKM-unitarity. On the other hand, it is worth noting 
that this value for brings the Big Bang nucleo-synthesis scenario in better agreement 
with the independent determination of the baryon content of the Universe deduced from 
the WMAP analysis of the cosmic microwave background power spectrum [28]. 

A new generation of r„-experiments, performed with magnetic storage devices where 
wall contacts of neutrons and thus neutron losses are avoided, is under construction. In 
this new approach, neutrons are trapped magnetically with permanent magnets [29] or 
with the superconducting magnets [30,31] (see also the proceedings of [12,25]). 

2.2.3 Results 

Inverting Eq. © leads to the following master formula [6,9]: 


,2 (4908.7 ±1.9)s 

\Vud\ — - 


( 10 ) 



Employing = 885.7(7)s and A = 1.2720(18) this implies 

\Vud\ = 0.9730 ± 0.0004 ± 0.0012 ± ±0.0002 , 


( 11 ) 


where the errors stem from the experimental uncertainty in the neutron lifetime, the /3- 
asymmetry Aq and the theoretical radiative corrections, respectively. As can be noted, so 
far the error is dominated by experimental uncertainties. A lower lifetime = 878.5(8)s 
would change \Vud\ to 0.9769(13). 

2.3 Vud from Tifj decays: the PIBETA experiment 

Pion beta decay, 7r+ —>• (also labeled tt^ and vTes), provides a theoretically excep¬ 

tionally clean means to study weak u-d quark mixing, i.e., the GKM matrix element Vud- 
Recent calculations of radiative corrections in Ref. [32] and Ref. [33] demonstrate that 
the theoretical uncertainty accompanying Vud extraction from pion beta decay is below 
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Figure 1; Schematic cross sec¬ 
tion of the PIBETA apparatus 
showing the main components: 
beam counters (BC, ACl, 
AC2), active degrader (AD) 
and target (AT), MWPC’s 
and support, plastic veto (PV) 
detectors and PMT’s, pure Csl 
calorimeter and PMT’s. 


0.05%. The new theoretical analysis of Ref. [9] indicates a further improvement in theo¬ 
retical accuracy. In the past, the low branching ratio of pion beta decay, rsj 10 has 
stood in the way of fully exploiting this opportunity experimentally. Until recently, the 
most accurate published measurement of was the one made by McFarlane et al. [34], 
with relative uncertainty AB/B barely under 4%, not even low enough to test the validity 
of radiative corrections which amount to approximately 3.2%. 

The PIBETA collaboration^ has initiated a program of measurements to improve the 
experimental precision of the branching ratios of pion beta and other rare pion as well as 
muon decays at the Paul Scherrer Institute (PSI) in Switzerland. Detailed information on 
the experimental method is available in Ref. [35] and at http://pibeta.phys.virginia.edu/. 

The experiment measures decays at rest. A 114MeV/c pion beam is tagged in an 
upstream beam counter (BC), slowed down in an active degrader (AD), and stopped in 
a segmented 9-element active target (AT). Charged particles are tracked in a pair of thin 
concentric multiwire proportional chambers (MWPC) and a 20-bar thin plastic scintillator 
veto hodoscope (PV). A large-acceptance 240-element spherical electromagnetic shower 
calorimeter made of pure Csl (12 radiation lengths thick) is used to detect the energy of 
both charged and neutral particles. Layout of the experiment is depicted schematically 
in Fig. n The 7r~^ —*■ decay events {7^e2) were used for normalization. 

The first phase of the experiment ran from 1999 to 2001. All observed decays were 
analyzed in parallel in order to minimize systematic uncertainties through multiple re¬ 
dundant crosschecks. A follow-up run, optimized for the tt’*' —radiative decay, took 
place in the summer of 2004. The main results of the 1999-2001 run were published in 


^Univ. of Virginia, PSI, Arizona State Univ. , JINR Dubna, Swierk, Tbilisi, IRB Zagreb 


























Refs. [36,37]. Because of the normalization to 7re2 decays, one can evaluate the pion beta 
decay branching ratio in two ways. Assuming full validity of the Standard Model, we can 
multiply the measured ratio i?7r/3/7?7re2 by the theoretical value of R™ = 1.2352 (5) x 10“^ 
from Ref. [38] and obtain (method 1): 

= [1.040 ± 0.004 (stat) ± 0.004 (syst)] x 10“® . (12) 

Alternatively, taking the experimental world average for the 7re2 branching ratio, = 
1.230 (4) X lO”'^ [22], one obtains (method 2): 

= [1.036 ± 0.004 (stat) ± 0.004 (syst) ± 0.003 (vTea)] x 10"® . (13) 

In either case, the result is in excellent agreement with the 90% C.L. SM limits: Bt^is = 
1.038 — 1.041 X 10“® [22]. These results represent the most stringent test of the CVC 
hypothesis in a meson performed so far. Furthermore, they lie about 6cr outside the SM 
limits if radiative corrections are excluded, thus conhrming the latter for the first time. 
The value of Vud obtained from these results: 14d = 0.9748 (25) (method 1), or 0.9728 (30) 
(method 2) is not yet competitive, but it agrees well with the one of superallowed transi¬ 
tions. We stress that this is not the hnal result of the PIBETA collaboration: an update 
with improved accuracy will be forthcoming. However, bringing the error down by an 
order of magnitude will require a whole new phase of the experiment. 

In the next stage of the project, the PIBETA collaboration will turn its attention to 
a new precise measurement of the branching ratio. 


3 Vus from K decays 


3.1 Theoretical aspects of decays 


The decay rates for all four modes {K = K^, i = fi,e) can be written compactly 

as follows: 


r(^o[7i) = 128^3 


x/f"-(0)|^x 


1 + 2 ^SU{2) 


2 

^ ^EM 


(14) 


Here Gp is the Fermi constant as extracted from muon decay, Aew = 1 + ^ (l — x 
log ^ -|- O(y^) represents the short distance electroweak correction to semileptonic 

charged-current processes, is a Clebsh-Gordan coefficient equal to 1 (1/a/2) for neu¬ 
tral (charged) kaon decay, while is a phase-space integral depending on slope and 

curvature of the form factors. The latter are defined by the QCD matrix elements 

{'K\Pn)\si^u\K\pK)) = (t) (Pi^ +P^)m + {PK-Pn)t, ■ (15) 


In the physical region these can be conveniently parameterised as 


Jo 


it) 


-/+,0 


(*) 






/+ " ( 0 ) ( 1 + + • • • ) > 


Mi 


(16) 

(17) 
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where t = {px — PttY- 

As shown explicitly in Eq. (d, it is convenient to normalise the form factors of all 
channels to (0), which in the following will simply be denoted by /+(0). The channel- 
dependent terms A|p^ 2 ) represent the isospin-breaking and long-distance elec¬ 

tromagnetic corrections, respectively. A determination of V^s from Ki^ decays at the 1% 
level requires ~ 1% theoretical control on /+(0) as well as the inclusion of A|p( 2 ) ^ind 

^EM- 

3.1.1 SU(2) breaking and radiative corrections 

The natural framework to analyze these corrections is provided by chiral perturbation 
theory [39-41] (CHPT), the low energy effective theory of QCD. Physical amplitudes 
are systematically expanded in powers of external momenta of pseudo-Goldstone bosons 
(tt, K, p) and quark masses. When including electromagnetic corrections, the power count¬ 
ing is in (e^)™ (p^/A^)"', with A^^ ~ and p^ ~ 0{pl^^, ~ 0{mg). To a given 

order in the above expansion, the effective theory contains a number of low energy cou¬ 
plings (LECs) unconstrained by symmetry alone. In lowest orders one retains model- 
independent predictive power, as these couplings can be determined by htting a subset of 
available observables. Even in higher orders the effective theory framework proves very 
useful, as it allows one to bound unknown or neglected terms via power counting and 
dimensional analysis arguments. 

Strong isospin breaking effects 0{niu — were hrst studied to O(p^) in Ref. [42]. 
Both loop and LECs contributions appear to this order. Using updated input on quark 
masses and the relevant LECs, the results quoted in Table El for A|p^ 2 ) were obtained in 
Ref. [43]. 

Long distance electromagnetic corrections were studied within CHPT to order 
in Refs. [43,44]. To this order, both virtual and real photon corrections contribute to 
Afi^. The virtual photon corrections involve (known) loops and tree level diagrams 
with insertion of O(e^p^) LECs. Some of these LECs have been estimated in [45], using 
large-A'c' techniques. The remaining LECs have been bounded by dimensional analysis 
(although they could be in principle estimated with the same techniques as in [45]). The 
resulting uncertainty is reported in Table El and does not affect the extraction of Vus at 
an appreciable level. 

Radiation of real photons is also an important ingredient in the calculation of A^j^, 
because only the inclusive sum of and rates is infrared hnite to any order in 
a. Moreover, the correction factor depends on the precise dehnition of inclusive rate. 
In Table El we collect results for the fully inclusive rate (“full”) and for the “3-body” 
prescription, where only radiative events consistent with three-body kinematics are kept. 
CHPT power counting implies that to order e^p^ one has to treat K and tt as point- 
like (and with constant weak form factors) in the calculation of the radiative rate, while 
structure dependent effects enter to higher order in the chiral expansion [46]. 

Radiative corrections to decays have been recently calculated also outside the 
CHPT framework [47,48]. Within these schemes, the UV divergences of loops are regu- 


10 



^SU{2) (^) 

AS(%) 

3-body full 


2.31 ± 0.22 [42,43] 

-0.35 ± 0.16 [43] 

-0.10 ± 0.16 [43] 


0 

+0.30 ± 0.10 [44] 

+0.55 ± 0.10 [44] 

+0.65 ± 0.15 [47] 


2.31 ± 0.22 [42,43] 




0 


+0.95 ± 0.15 [47] 


Table 2: Summary of SU(2) and radiative correction factors for various Ki^ decay modes. 
Refs. [42-44] work within chiral perturbation theory to order e^p^, while Ref. [47] works 
within a hadronic model for Kaon electromagnetic interactions. 


lated with a cutoff (chosen to be around 1 GeV). In addition, the treatment of radiative 
decays includes part of the structure dependent effects, introduced by the use of form 
factors in the weak vertices. Table |21 shows that numerically the “model” approach of 
Ref. [47] agrees rather well with the effective theory. 

Finally, it is worth stressing that the consistency of the calculated strong and SU(2) 
corrections can be probed by experimental data by comparing the determination of the 
universal term Vus x /+(0) from the various decay modes, as we will discuss in section l!01 


3.1.2 Analytic results on /+(0) 

Within CHPT we can break up the form factor according to its expansion in quark masses: 


/+( 0 ) — 1 + /p4 + /p 6 + . . . . (18) 

Deviations from unity (the octet symmetry limit) are of second order in SU(3) breaking [4]. 
The hrst correction arises to O(p^) in CHPT: a hnite one-loop contribution [42,49] de¬ 
termines /p 4 = —0.0227 in terms of TV, Mk and M^r, with essentially no uncertainty. 
The p® term receives contributions from pure two-loop diagrams, one-loop diagrams with 
insertion of one vertex from the effective Lagrangian, and pure tree-level diagrams with 
two insertions from the Lagrangian or one insertion from the p® Lagrangian [50,51]: 

+ /•;“(;.). ( 19 ) 

Individual components depend on the chiral renormalization scale p, their sum being 
scale independent. Using as reference scale p = Mp = 0.77 GeV and the Li from £t 10 in 
Ref. [53], one has [51]: 

fyoops^Mp) = 0.0113 , = -0.0020 ± 0.0005 . (20) 

The explicit form for the tree-level contribution in terms of the LECs L 5 [41] and ^ 12,34 [54] 
is then [51,52] 


fXiMp) = 8 


(Ml - Ml 
F3 


F3 


- Cl,{Mp) - c;,{Mp) 


( 21 ) 
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Ll{Mp) is known from phenomenology to a level that induces less than 1 % uncertainty 
in The p® constants 6 * 12,34 can in principle be determined phenomenologically. 

It has been shown in Ref. [51] that combinations of C 12 and 6*34 govern the slope Aq 
and curvature Ag of the scalar form factor /o(t), accessible in i 6^3 decays. In order to 
extract C 12 and 6*34 to a useful level (i.e. leading to 1 % hnal uncertainty in 14 s), one 
needs experimental errors at the level AAg ~ 0.001 (roughly a 5% measurement) and 
AAg ~ 0.0001 (roughly a 20 % measurement), as well as Fk/F^^ at the 1 % level from 
theory. 


harge-Nc estimate of 

In Ref. [55] a (truncated) laige-Nc estimate of was performed. It was based on 
matching a meromorphic approximation to the (SPP) Green function (with poles corre¬ 
sponding to the lowest-lying scalar and pseudoscalar resonances) onto QCD by imposing 
the correct large-momentum falloff, both off-shell and on one- and two-pion mass shells. 
In particular, G 12 is uniquely determined by requiring the correct behavior of the pion 
scalar form factor (7r|S'|7r), while G 34 is hxed by the correct scaling of the one-pion form 
factors {7i\S\V) and (7r|P|5). The uncertainty of the large-i\4 matching procedure was 
estimated by varying the chiral renormalization scale at which the matching is performed 
in the range fi G [M^, IGeV], and is found to be = ±0.008. The hnal result is 


fTiMp) = 


{Ml - Mir 
2 Ml 


Ml 


= — 0.002 ± 0.008 ± 0.002 


( 22 ) 


and is much smaller than the ratio of mass scales {Ml — M^Y/Ml would suggest, due to 
interfering contributions. When combined with the p® loop corrections [51], this estimate 
leads to fp6 = 0.007 ± 0.012. Variations of the hadronic ansatz lead to the conclusion 
that the smallness of the tree-level part compared to the loop contribution of 0 (p®) for 
/+( 0 ) appears as generic feature of a few-resonance approximation for the set of large- 
momentum constraints considered. As a consitency check of this approach, it is worth 
to mention that within the same framework one obtains a prediction for the slope of the 
scalar form factor, Aq = 0.013 ± 0.002 i/^vc ± 0.001 fully consistent with the value 
measured by KTeV in decays, Aq = (13.72 ± 1.31) x 10“^ [56]. 

Combining this p® result with the well-know p^ term, leads to the following global 
estimate of /+( 0 ): 

/+(0)large-iVe = 0.984 ± 0.012 . (23) 

This value is substantially higher -although compatible within the errors- with respect 
to the old estimate by Leutwyler Roos [49] 


/+(0)Leutwyler-Roos = 0.961 ± 0.008 , (24) 

which for a long time has been the reference value of /+( 0 ) (and it is still the value adopted 
by the PDG [22]) in the extraction of 14s- 
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3.1.3 /+(0) from lattice QCD 

Starting from Ref. [57], in the last few months it has been realized that lattice QCD is 
an excellent tool to estimate /+(0) at a level of accuracy interesting for phenomenological 
purposes [58-62], 

On general grounds, determining a form factor at the 1% level of accuracy seems 
very challenging -if not impossible- for present lattice-QCD calculations. However, the 
specihc case of /+(0) is quite special: by appropriate ratios of correlation functions one 
can directly isolate the SU(3)-breaking quantity [/+(0) — 1], or even better the quantity 
[/+(0) — 1 — fp4,], which is the only irreducible source of uncertainty [57]. Estimating these 
SU(3)-breaking quantities with a relative error of about 30% is sufficient to predict /+(0) 
at the 1% level or below. Thus even with the present techniques there are good prospects 
to obtain lattice estimates of /+(0) of phenomenological interest. 

The analysis of Ref. [57] is based on the following three main steps: 


I Evaluation of the scalar form factor fo{q‘^) at = {Mk — Applying a 

method originally proposed in Ref. [58] to investigate heavy-light form factors, the 
scalar form factor is extracted from the following double ratio of matrix elements: 


{7i\s^ou\K){K\u^os\7i) 

{7i\u'you\7i){K\s'yos\K) 


{Mk + 
4:MkM^ 




(25) 


where all mesons are at rest. The double ratio and the kinematical configuration al¬ 
low to reduce most of the systematic uncertainties and to reach a statistical accuracy 
on /o(Q'max) below 1% (see figure |2I). 

2 Extrapolation of /o(Q'max) /o(0) = /-i-(0). By evaluating the q^ dependence of 
the form factor, the latter is extrapolated from q^^.^ to q^ = 0. This procedure 
is performed independently for various sets of meson masses (with corresponding 
light-quark masses chosen in the range Q.dm^- 2ms), and using different functional 
forms (linear, quadratic and polar) for the dependence. A byproduct of this step 
is an estimate of the physical slope of the scalar form factor: the prediction quoted 
in Ref. [57] for this observable has recently been confirmed by experimental analysis 
of KTeV [56]. 


3 Subtraction of the chiral logs and chiral extrapolation. The /+(0)-values thus ob¬ 
tained needs to be extrapolated to the physical values of Mk and M^,. In order to 
reduce the error of this extrapolation, and correct for the leading quenched artifacts, 
the following ratio is considered: 


R{Mk,M^) 


A/ 

(AM2)2 


1 + /;.(A'/s-, - /+( 0 ; Mk, M„) 

(Ml - Mi? 


(26) 


Here f^^{MK, M^,) denotes the 0{p'^) contribution evaluated within quenched CHPT 
which, similarly to its unquenched analog {fp^), is finite and free from unknown 
counterterms. By construction, the ratio (|2ni) is finite in the SU(3) limit, does not 
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Figure 2: Left: /o(5'max) ^ function of [M]^ — M^] in units of the lattice spacing (a). 

Right: R{Mk, as a function of [Mj^ + M^] for the cases of the linear, quadratic and 
polar g^-extrapolation in the £t of /o(q'^); the triangles indicate the values of R{Mk, 
extrapolated to the physical kaon and pion masses. 


depend on any subtraction scale, and is free from the dominant quenched chiral 
logs of 0{p'^). Extrapolating R{Mk, M^-) to the physical masses (see hgureH)) leads 
to A/ = X [(AM2)2]phy« = (0.017 ± O.OOSstat ± O-OOT^yst), which 

implies [57]: 

/+(O)Lattice—quenched 0.960 ± O.OOSstat ± 0.007,yst , (27) 

in remarkable agreement with the old estimate by Leutwyler and Roos [see Eq. ()24|l ]. 

The result of Ref. [57] is affected by an irreducible systematic error due to the residual 
effects [0{p^) and above] of the quenched approximation. This error is difficult to quan¬ 
tify and is not included in the systematic error of Eq. dZIl), which is dominated by the 
extrapolation to the physical meson masses. However, the good agreement between the 
slope of the scalar form factor obtained within this simulation and the experimental one 
seems to suggest that this additional source of uncertainty is not large. 

Very recently new analyses of /+(0) with unquenched lattice simulations {Nf = 2 and 
Nf = 2-1-1) have been started by various groups [59]. In particular, the preliminary 
results by the collaborations MILC and HPQCD (/+(0) = 0.962(6)(9) [60]), JLQCD 
(/+(0) = 0.952(6) [61]), and RBC (/+(0) = 0.955(12) [62]) are very encouraging: the 
errors do not exceed the 1% level and the results are well compatible. Interestingly, these 
preliminary unquenched results for /+(0) are also in good agreement with the quenched 
estimate of Ref. [57]. 
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3.2 The extraction of V^s from decays 

In the last few years several experiments, based on different techniques, have put a lot 
of effort in re-measuring all the observables of the main K decay modes, whose values 
available in literature dated back to the 70’s. As discussed above, the master formula 
for the extraction of |I4s| from Ki^ decays is Eq. fll4D . The experimental inputs are the 
semileptonic widths (based on the semileptonic branching fractions and lifetimes) and the 
form factors, which are necessary for the calculation of the phase-space integrals. Since 
none of the experiments have measured yet all of the experimental inputs required to 
calculate |14<j| independently, we have calculated average values of |I4s|/-i-(0) for A'+(e3), 
Kl^cS), KL^fiS), and Ks{e3) using the inputs described in the following sections. We 
then combine these results to hnd an average |I 4 s|/+( 0 ) for decays. 

3.2.1 Charged and neutral kaon branching fractions 

Three experiments have contributed to the measurement of the neutral kaon branching 
fractions: KTeV, NA48 and KLOE. The hrst two are hxed target detectors using sec¬ 
ondary high energy neutral beams and situated downstream from the primary target to 
obtain a beam of Kl, with small contamination of Ks- No absolute kaon count is available 
so that one can measure ratios of partial widths. Instead KLOE is a collider experiment 
at a 0 -factory where Ks — Kl pairs are produced and therefore takes advantage of the 
unique feature of the tagging: identihed Kl{Ks) decays tag a Ks{Kl) beam and provide 
the means for measurements of absolute branching ratios. Therefore the measurements 
performed by the three experiments have quite different systematic uncertainties. 

For the Kl branching fractions, we consider the following experimental inputs: 

• KTeV measured the following 5 partial width ratios [65]: 

T{Kl —> gKv)/T{K l —> , T{Kl —> 7r+7r“7r°)/r(iLi —>• vr^e^z/), 

T{Kl —> /T{Kl —>■ TT^e^z/) , T{Kl —>■ /T{Kl —>■ vr^e^z/) , 

V[Kl —> /V[Kl —> 

Since the six decay modes listed above account for more than 99.9% of the total 
decay rate, the hve partial width ratios may be converted into measurements of the 
branching fractions for the six decay modes. 

• KLOE uses a tagged Kl sample to measure the 4 largest Kl branching fractions [ 66 ]. 

• NA48 measures the following 2 ratios [ 68 ]: 

T{Kl —> 'K^e^v)/V{KL —> 2 tracks), T{Kl —>■ 7r°7r°7r°)/r(JL5 —> tt^tt^), 
which can be used to determine B{Ke^) and i?( 37 r°). 

A critical issue is represented by the inclusiveness of the measurement and the treatment 
of the radiative photon. These have been carefully accounted for in all the experiments. 
A £t to all of these measurements, accounting for correlations, gives the Kl branching 
ratios summarized in Table El Figure El shows a comparison of these measurements along 
with the best £t values for each of the six branching fractions. 
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Decay Mode 

Branching fraction 

Pi (lO^s-^) 

Kl- 

->■ TT^e^Z/ 

0.4040 ± 0.0008 

0.7908 ±0.0032 

Kl- 

H- 

-H 

t 

0.2699 ± 0.0008 

0.5283 ±0.0023 

Kl- 

->■ TT+TT 7r° 

0.1253 ±0.0006 

0.2452 ±0.0015 

Kl- 

0 

0 

0 

t 

0.1972 ±0.0012 

0.3859 ±0.0029 

Kl- 

1 

+ 

t 

(1.971 ±0.012) X 10-3 

(3.857 ±0.027) x 10-3 

Kl- 

0 

0 

t 

(0.880 ±0.008) X 10-3 

(1.722 ±0.017) X 10-3 


Table 3: Average Kl branching fractions and widths based on £t to all new measnre- 
ments from KTeV, KLOE, and NA48. The width measurements use the PDG 2004 
average lifetime combined with the two preliminary KLOE measurements mentioned 
in Section 13.2.21 

For Ks, we use the KLOE preliminary measurement: B{Ks —> vrez/) = (7.09 ±0.09) x 
10“^ [71]. For we average the BNL E865 [72] measurement, B{K^ —> vr^e^z/) = 
(5.13 ± 0.10)%, with the more recent NA48 result, B{K^ —> 7 r°e^z/) = (5.14 ± 0.06)%, 
Ending B{K^ —> ■n^e^v) = (5.14 ± 0.06)%. 

3.2.2 Lifetime Measurements 

KLOE presented two new measurements of the lifetime: an “indirect method,” based 
on the Kl lifetime required to make = 1 , exploiting the lifetime dependence of the 
detector acceptance [ 66 ], and a “direct method”, based on the fit of the proper time 
distribution of JLi, —> 37 r° decays [73]. These new results and the old PDG average are 
listed in Table El The new average value, which we use for the results quoted below, is 
Ti = (50.98 ±0.21) ns. 

Gombining the Ki branching fractions with the new lifetime gives the partial decay 
widths quoted in Table [31 Note that correlations between the KLOE branching fractions 
and the “indirect” KLOE lifetime determination have been included. For the Ks and K^ 
lifetimes, we use the PDG averages. 

It is worthwhile to note that KLOE is about to report new results on the K^ lifetimes 
and the achievable accuracy is at the few per mil level. The measurement is done with 
two independent methods: the first one is based on the charged kaon decay path while 
the second one measures directly the kaon time of flight using the photons from vr® — 77 
decay, using the charged kaon decay channels with a vr^ in the final state. This will allow 
the clarification of the present situation which shows discrepancies between “in-flight” 
and “at-rest” charged kaon lifetime measurements, used in the PDG average. 

3.2.3 Phase Space Integrals 

Recent experiments have also performed new measurements of the semileptonic form 
factors needed to calculate the phase space integrals: 

• KTeV has measured form factors in both K^^ and [67]. 


16 









N_A48 






KIEV 


0.38 0.39 0.4 

B(KL->7iev) 



KLOE 



f<TEV 
1-•— 

'pDG 02 




0.41 


0.27 

B(Kl->ji|j.v) 


N 

^48 


A 






0.19 


0.2 


0.21 


0.124 


0 . 126 „ 


B(Kl^3ji^) 




KIEV 
1 


PDG 02 
I-O-1 


1.9 2 2.1 

1000xB(K,^jiV) 


0.275 


1 — 


:loe 

K+EV 
—•- 1 

PDG 

02 





0.128 


KTE 

/ 




PDG 02 

1-0-1 


0.85 0.9 0.95 

1000X B(Kl^jiV) 


Figure 3: Kl branching fractions measured by KTeV, KLOE, NA48, and from the PDG 
2002 £t (open circles). The vertical lines indicate the ilcr bounds from a fit to all KTeV, 
KLOE, and NA48 measurements. 

• NA48 has measured the form factor [69]. 

• ISTRA+ has measured the K~^ form factor [70]. 

KLOE is also measuring the form factors in both neutral and charged kaon decays, al¬ 
though final results have not been published yet. Note that, in principle, KLOE is the 
only experiment with the possibility of measuring all the useful inputs for the extraction 
of Vus, namely lifetimes, branching fractions and form factors both for neutral and charged 
kaons. 

In the present analysis, to calculate phase space integrals, we use the KTeV quadratic 
form factor results for neutral kaon decays and the ISTRA-I- quadratic form factor mea¬ 
surements for charged kaons. For both charged and neutral decays, we include an ad¬ 
ditional 0.7% uncertainty to the phase space integrals, as suggested by KTeV [67], to 
account for differences between the quadratic and pole model form factor parameteriza- 
tions, both of which give acceptable fits to the data. The resulting phase space integrals 
are = 0.1535 ± 0.0011, = 0.10165 ± 0.0008, and = 0.1591 ± 0.00012. 

3.2.4 JG 3 results for |14s| 

We are now ready to extract four estimates of /+(0)|14s| from Eq. (fTlj) . using the available 
data on Kl (both e and /i modes), and Ks decays. To this purpose, we set Sew = 
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Source 

Lifetime (ns) 

PDG 2004 Average 

51.5 ±0.4 

KLOE (“indirect”) 

50.72 ±0.35 

KLOE (“direct”) 

50.87 ±0.31 

New Average 

50.98 ±0.21 


Table 4: Kl lifetime measurements. 


1.023 [74] and use the SU(2) and radiative correction factors reported in Table El 

Figure Elshows a comparison of the PDG and the averages of recent measurements for 
f+{0)\Vus\ for the various decay modes. The hgure also shows /+(0)(1 — |14dP — 
namely the expectation for /+(0)|14s| assuming unitarity, based on \ Vud\ = 0.9738±0.0003, 
iKifel = (3.6 ± 0.7) X 10“^, and several recent calculations of /+(0). The average of all 
recent measurements gives 


f+{0)\Vus\ = 0.2173 ± 0.0008 (28) 

Using the Leutwyler-Roos estimate of /+(0), this implies 

= 0.2261 ±0.0021 [/+(0) = 0.961 ± 0.008] . (29) 


3.3 Vus from K ^2 

As recently pointed out in Ref. [76], besides the traditional method of extracting |14s| 
from Ki 2 , decays, one can obtain an independent and competitive estimate of | 14 s| (or, to 
be more precise, of | 14 s/hud|) from the ratio of the inclusive decay widths of > /i+i/( 7 ) 
and decays. The measurable ratio can be written as follows 


r{K* ^ 7, Fj Ml(Ml - Mlf r g _ 

r(7r+ ^ ,^+>.(7)) VS Fi MJ(MJ -MlfV 7 ^ 


(30) 


where Fk^-k denote kaon and pion decay constants, and Ck,-k parametrize the radiative- 
inclusive electroweak corrections, accounting for soft-photon emission and corresponding 
virtual corrections. According to the detailed analyses of Refs. [78,79]: — Ck = 

3.0 ±0.75. 

The key theoretical input for this method is the estimate of Fk/F.,^, obtained by means 
of lattice QCD. Contrary to the case of /+(0), which is protected by the Ademollo-Gatto 
theorem, the quantity [Fk/F^^ — 1] breaks SU(3) invariance already at the hrst order. It 
is therefore very challenging to estimate F^/F^^ at the 1% level of accuracy. However, 
the MILC collaboration has recently shown that this is possible [80]. In 2004 they found 
Fk/Ft^ = 1.210 ±0.014 [80], value which is compatible with the preliminary result of their 
new analysis, namely FV/FV = 1.198 ± 0.003](o;oo5 [81]. 

On the experimental side, the dominant error in Eq. iHD is induced by the kaon decay 
width. KLOE has recently performed a new measurement of the corresponding absolute 


18 












“1-1-1-1-1-1-1-1-1-1-r 

I-B- 1 


T-1-1-1-1-r 

—B-1 


1-1-1-1-1-r 


K 


+ 


K"e3 

[PDG 2002] 

K"e3 

[new] 

KLe3 

[PDG 2002] 

KLe3 

[new] 

Kl,3 

[new] 

Kse3 

[new] 


I-B-1 


K, 




f,(0) (l-IV,J2-IV„/)‘'- 

Leutwyler Roos 
Lattice [Becirevic et al.] 
Large [Cirigliano et al.] 


_l_I_I_I_L_ 




_l_I_I_I_ I _L_ 


-e-1 


e- 


h 


J_ 


-e- 


J_I_I_I_ I _I_I_I_L_ 


0.205 


0.21 


0.215 0.22 

f,(0) lv,J 


0.225 


0.23 


Figure 4: Comparison of the recent measurements of \Vus\f+{0) with the PDG 2002 
results, and with the expectations of unitarity based on different theoretical calculations 
of /+(0) [49,55,57]. The values of /+(0)(1 — |14dP — have been obtained using 

the 2005 value of Vud in Eq. the error bars are largely dominated by the theoretical 
uncertainty on /+(0). The vertical lines correpsond to the average value of /+(0)|14s| in 
Eq. (j2Hl). 

branching ratio, obtaining B{K^ /i+i/( 7 )) = 0.6366 ± 0.009 ± 0.00015 [77]. This mea¬ 
surement is fully inclusive of the hnal state radiation and is based on a sample of 
events tagged by K~ decays produced at the 0-meson peak, thus avoiding normaliza¬ 
tion issues such as trigger and reconstruction efficiency (which enter at hrst order in the 
determination of the branching ratio). 

Using the new KLOE result, the life time from PDG and Vud = 0.9738 ± 0.0003 
(see section ini), leads to 

\Vus\K+^^+.i^) = 0.2223 ±0.0026 = 1.210 ± 0.014) , (31) 

= 0.2245l‘];°[]^} = 1.198 ± 0.003^^© , (32) 

In both cases the error is largely dominated by the uncertainty on /F^- and, in partic¬ 
ular, the uncertainty induced by 110^1 is negligible. Both results are consistent with the 
Vus determination from 
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4 Alternative approaches to Vus 


4.1 Vus from hyperon decays: status and perspectives 

Hyperon semileptonic decays (HSD), denoted here as Bi —> B 2 + V + Vi^ are described by 
the V — A theory. This theory states that the weak transitions arise from the self-coupling 
of a single charged current, which is the sum of the current operators for the leptons and 
the strongly interacting particles. Each current can be expressed as a linear combination 
of vector and axial terms. The matrix elements of the lepton current are unambiguous. 
In contrast, non-perturbative strong interaction effects at low energy weak interactions 
force the introduction of phenomenological form factors to account for strong interaction 
effects in the matrix elements. The transition amplitude for HSD can be written as [82] 

Mo = - 15 )Vu{Pu)][uB2{P2)WuUb,{Pi)], (33) 


where Vckm is either Vud or 14s and 


Wu = fi{q^hu + 


Ml 




fsjq^) 

Ml 


qfi + 


9liq‘^h^. + 


92{q^ 

Ml 


-(y^luq'' + 


9^{q^ 

Ml 




75- (34) 


Here g = pi — p 2 is the momentum transfer and Mi is the mass of the decaying hyperon. 
The quantities fi{q^) and gi{q^) are the conventional vector and axial-vector form factors, 
which are required to be real by time reversal invariance. 

The usual approach to evaluate the expected properties of the weak hadronic current 
is based on the flavor SU(3) symmetry of the strong interactions. In the limit of exact 
SU(3) symmetry, the hadronic weak currents belong to SU(3) octets, so that the form 
factors of different HSD are related to each other. The vector part of the weak current 
and the electromagnetic current belong to the same octet. Thus, the vector form factors 
are related at = 0 to the electric charge and the anomalous magnetic moments of the 
nucleons. The conservation of the electromagnetic current implies that fo{q^) vanishes 
for all HSD in the SU(3) limit. Similarly, gi is given in terms of two reduced form 
factors F and D whereas g 2 for diagonal matrix elements of hermitian currents vanishes 
by hermiticity and time-reversal invariance. SU(3) symmetry then implies that ^2 = 0 in 
the symmetry limit. 

Currently, experiments in HSD allow precise measurements of form factors (for a re¬ 
view about the experimental situation of HSD see for instance Ref. [83]). The statistical 
errors of these experiments are rather small, and more effort has been put into the reduc¬ 
tion of the systematic errors, which can be of two types. The hrst one comes from the 
different shortcomings of the experimental devices. The second one, of theoretical nature, 
comprises i) radiative corrections and ii) theoretical assumptions for some form factors, 
including their g^ dependence. Let us briefly discuss each one. 

The analysis of low, medium, and high-statistics experiments of HSD requires the in¬ 
clusion of radiative corrections. Since no hrst principle calculation of radiative corrections 
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is yet possible, these corrections are committed to model dependence and the experimen¬ 
tal analyses which use them become mo del-dependent too. Up to order (a/7r)(g/Mi), the 
model dependence of radiative corrections can be absorbed into the form factors origi¬ 
nally dehned in the matrix elements of the hadronic current [82] and the remaining 
part is model independent. Within these orders of approximation one is left with general 
expressions which can be used in model-independent analyses [84]. Among the integrated 
observables in HSD, only the decay rates need to be corrected whereas the angular corre¬ 
lations and spin-asymmetry coefficients are practically unaffected [82]. 

On the other hand, the assumptions on the form factors are subtle to handle. In 
particular, their g^-dependence cannot always be neglected, since noticeable contributions 
can arise. In order to obtain expressions correct to order (9(g^), the g^-dependence of /2 
and g 2 can be ignored, because they already contribute to order 0{q) to the decay rate. 
For /i(g^) and gi(g^), instead, a linear expansion in g^ must be considered, because higher 
powers amount to negligible contributions to the decay rate. A dipole parametrization of 
the leading form factors works hne [82]. 

Another important issue that must be taken into account is the validity of the exact 
SU(3) limit. Presently, experiments are precise enough that exact SU(3) no longer yields 
a reliable £t. For HSD, due to the presence of the axial current, SU(3) breaking (SB) 
can occur in hrst order. This fact makes HSD to be apparently less reliable to use for 
determining Vus than Ki^ decays, where the Ademollo-Gatto theorem reduces the effects 
of SB. For this reason the present value of Vus quoted by the Particle Data Group [22] is 
essentially derived from Ki^, while the one from HSD is discarded. 

Although there are various treatments on the calculation of SB in HSD [85-89], it is 
hard to assess their success because their predictions vary substantially from one another 
and none of them can be considered as fully consistent. These calculations incorporate 
second-order SB corrections into /i. Some computations [85-87] hnd that /i, for | AAI = 1 
decays, is reduced from its symmetry limit value such that < 1, whereas 

the calculation of Ref. [88], performed in the framework of Heavy Baryon Chiral Per¬ 
turbation Theory (HBChPT) [90], found large and positive corrections. Recently, the 
calculation of SB corrections within HBChPT has been reconsidered in Ref. [89], where 
it is found that the calculation of [88] contained some mistakes and missed some im¬ 
portant contributions. The result of Ref. [89] is that SB corrections, as obtained from 
HBChPT, are still positive but smaller than previously claimed in [88]. In addition, the 
quantitative estimate of these corrections strongly depends on the values used for the low 
energy constants and does not include the contribution coming from the baryon decup let. 
This suggests that the corrections coming from higher orders can be large, and that the 
present estimates may not be reliable. Alternatively, SB effects can be also extracted 
from the data, for instance parameterising them in the framework of the l/Nn expansion 
of QCD [91-94]. 

Fits to the experimental data of HSD to extract 14s can be performed using the 
decay rates and the spin and angular correlation coefficients [22]. There are sufficient 
data from hve decays to make it possible: A —>• pe~Ve, ne~Ve, 2“ —>• Ae^W, 

S°e“z7e) and > S+e“z7e- An alternative set of data is constituted by the decay 
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rates and the measured values of the gi/ fi ratios. Since the latter, however, contain less 
experimental information, using the several angular coefficients instead of gi/ fi provides 
a more sensitive test. 

In order to extract Vus several analyses can be performed under different assumptions. 
The analysis done in Ref. [83] neglects the quadratic SB corrections in the vector form 
factor fi and accounts for the larger effects in the axial form factor gi by using the 
measured values of the gi/fi ratios of the above processes, except for E~ —>■ E°e“I7e- 
The matrix element Vus is then extracted separately in each decay and the results are 
combined to obtain the value Vus = 0.2250 ± 0.0027, which is in good agreement with the 
unitarity requirement. 

A similar analysis has been done in Ref. [94], this time using the decay rates and 
the angular correlation coefficients and performing a global ht to the data of the hve 
decays. By hrst assuming the validity of exact SU(3) symmetry, the analysis yields 
Vus = 0.2238 ± 0.0019, with a x^/dof of around 2.5. This high value of signal 

the presence of not negligible SB corrections. One can proceed further and study the 
effects upon Vus of the leading form factors when SB corrections are taken into account. 
If one hxes fi at their SU(3) symmetry values and incorporates hrst-order SB in gi one 
obtains Vus = 0.2230 ± 0.0019. Next, incorporating SB effects in both fi and gi yields 
Vus = 0.2199 ±0.0026 [94]. This analysis also hnds that SB corrections to fi increase their 
magnitudes over their SU(3) symmetric predictions by up to 7% and that corrections to 
gi are consistent with expectations. This latter Vus, although in good agreement with 
the one quoted by the Particle Data Group, does not help in improving the unitarity 
test. Another interesting hnding is that the incorporation of more rehned SB corrections 
systematically reduces the value of Vus from its SU(3) symmetric prediction, rather than 
increasing it to better satisfy unitarity. Of course an increase of Vus can be achieved with 
a symmetry breaking pattern such that < 1 [85-87], but current data do not 

seem to favour this trend [94]. 

We are now in a position to make some statements. The hrst one is that the assumption 
of exact SU(3) symmetry to compare theory and experiment in BSD is questionable due 
to the poor hts it produces. The second one is the fact that deviations from the exact 
SU(3) limit, in particular for the case of the leading axial form factor gi, are indeed 
important to reliably determine Vus from BSD, which can rival in precision with the one 
from Ki 3 decays. Nevertheless more work, theoretical and experimental, is required in 
the next future to be conclusive about Vus from BSD. Particularly, recent techniques in 
lattice QCD [95] allow to measure form factors with great accuracy, thus determining SB 
corrections in a model-independent manner. 

4.2 Determination of |Vus| and mg from hadronic r decays 

Already over a decade ago, it was realized [97] that hadronic r decays provide a very clean 
testing ground for low energy QCD. For example, analysing the non strange r spectral 
function [96] leads to a determination of the QCD coupling as that is competitive with 
the world average. Also, in view of the increasingly precise data on the strange spectral 
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function from ALEPH [98], OPAL [99] and CLEO [100], studying flavor breaking observ¬ 
ables allows for a determination of the strange quark mass, another fundamental QCD 
parameter. The first analyses of this kind [101-104] were plagued with large uncertainties 
in the scalar contributions due to a bad behavior of the associated perturbation series. 
This problem was circumvented [105] by substituting the theoretical expressions with phe¬ 
nomenological counterparts, thereby reducing the theoretical uncertainties significantly. 
Along with this procedure, it was realized that the main uncertainty in the strange quark 
mass determination was now arising from the uncertainty in the CKM parameter |Pus|, 
which suggests [105,106] to turn the analysis around and determine |14s|- 

The basic objects for a QCD analysis of hadronic r decays are the two point corre¬ 
lation functions for vector = qf^^qj and axial-vector AC = qf^^'^^qj currents. Both 
correlation functions contain (pseudo)scalar and (axial) vector contributions. Then, the 
ratio of hadronic to leptonic r decay rates 


Rn- = 


r[r- 


hadrons( 7 )] 


P [r ^ e ’ 

as well as higher moments of the invariant mass distribution 
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can be calculated in an operator product expansion: 
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Here, Sew = 1.0201 ± 0.0003 summarizes electroweak radiative corrections, while explicit 
expressions for the and additional theoretical information can be found in [101] 

as well as references therein. The dominant contribution to Rr is purely perturbative 
and summarized in while higher dimensional contributions that also depend on 

the flavor content are suppressed. The most important of these suppressed terms are 
proportional to and ms{qq)- To reduce perturbative uncertainties one conveniently 
discusses the flavor breaking observable 
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where the dehnition of the flavour-dependent moments and R^^'g can be found 

in Ref. [105]. In this expression, the main parametric uncertainties arise only from mg 
and |Ks|, so that one could ideally determine both parameters simultaneously from the 
experimental analysis of several moments (see [105] for more details). In absence of such 
an analysis, we proceed as in [106] and begin by determining |Ks| from the moment with 
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Moment 

(2,0) 

(3,0) 

(4,0) 

ms{Mr) MeV 

93.2+^^ 


79.2+^^ 


Table 5: Results for ms(M^) extracted from the different moments. 


the smallest dependence on i.e. the (0,0) moment. Using ms(2GeV) = (95 ±20) MeV 
(in the the MS scheme), in agreement with recent sum rule and lattice calculations, one 
hnds = 0.218 ± 0.026, and, with 


one obtains 


IK. 


2 
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( 0 , 0 ) 

T,S 


d(0,0) 

M,v+a 


-6R 


(0,0) 

T,th 


(38) 


IK.I = 0.2208 ± 0.0033exp ± 0.0009th = 0.2208 ± 0.0034, (39) 

where we have used Rfy]^j^ = 3.469 ± 0.014 and R'fs'^ = 0.167 ± 0.0050 from [99]. 
The most important feature of this determination is the small theoretical uncertainty, the 
reason for which can nicely be seen from Eq. (EHD: the large cancellations between strange 
and non strange channels lead to a small value of , so that the main sensitivity is 

to the experimental input. It is also interesting to check whether this value of |K.| is 
compatible with unitarity: using the value |Kd| = 0.9739 ± 0.0003 (see section 2), one 
hnds that Eq. m is violated only at the 1.8 a level. 

In the next step, one uses the value of |Ks| thus obtained and determines the strange 
quark mass from higher moments. One hnds the values given in Table El where the 
individual sources of the uncertainties are given in [106]. 

The weighted average of the strange mass values obtained for the diherent moments 
give 

ms{M^) = 84 ± 23 MeV ^ m,(2GeV) = 81 ± 22 MeV , (40) 

in good agreement with the average given above, while a more detailed comparison can 
be found in [106]. The values of rris display a monotonous k dependence, which has 
been signihcantly reduced with the new OPAL data, mainly due to the new value for the 
B{t~ —>■ branching ratio. 

In summary, the |Ks| value obtained from r decays is beginning to become competitive 
with the standard determination from K decays and is given in Eq. (P|) . In particular, 
due to the fact that the largest amount of the uncertainty is still experimental, one can 
expect that the better data samples from BaBar and Belle will reduce this uncertainty 
signihcantly. In view of these perspectives it will be important to have the possibility 
to determine both ms and |Ks| simultaneously. An analysis of this type is underway. 
Additionally, further, precise data will clarify whether the moment dependence of the ms 
results is indeed a purely experimental issue, and thereby allow a consistency check of the 
whole analysis. 
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Future prospects on Vcs and Vcd 

5.1 Theoretical developments 

As far as the extraction of Vcs and Vcd is concerned, the most interesting recent devel¬ 
opments are the lattice QCD calculations of the matrix elements for semileptonic and 
leptonic decays. With the advent of dynamical simulations with light sea quarks, 
decays to K and vr mesons are favoured because the hnal-state particles are stable to 
strong decay. The headline results have come from simulations using 2-1-1 flavours of im¬ 
proved staggered light quarks where the common light (m, d) quark mass is in the range 
m^/S < miight < 3ms/4. The “fourth root trick” used in generating these dynamical 
conhgurations has not been completely justihed theoretically, but on the other hand it 
has not so far failed any test. There is a requirement to deal with the unphysical “tastes” 
introduced, but this has been addressed by the development of staggered chiral pertur¬ 
bation theory (SyPT) [107,108], allowing the subtraction of discretisation effects arising 
from light quark taste violations. 

5.1.1 Lattice Results for Semileptonic D Decays 

For semileptonic D decays to a light hnal state pseudoscalar or vector meson, P or V, the 
squared momentum transfer is = m|, -|- mpy — 2m qE where E is the P or V energy in 
the D rest frame. In lattice simulations, the entire range, 0 < —mpyY = g^^x; 

can be accessed while keeping the initial and hnal state meson spatial momenta small 
enough to avoid large discretisation effects. 

Two form factors, /+(g^) and /o(g^) are needed to describe the vector current matrix 
element for semileptonic decays to a light pseudoscalar, but only /+(g^) is needed for 
the decay rate if the lepton mass can be ignored. Lattice QCD gives results for both 
form factors and this can be helpful since they can be ht simultaneously imposing the 
constraint /+(0) = /o(0). 

The Fermilab-MILC-HPQCD collaboration has reported results for semileptonic D 
K and D —»■ vr decays using 2-1-1 havours of improved staggered quarks [109,110]. The 
charm quark is implemented using the Fermilab method [111], while the calculations are 
so far performed at one lattice spacing. Indeed, the dominant systematic effect comes 
from the heavy quark discretisation error, estimated from the mismatch of the continuum 
and lattice heavy quark effective theories [112,113] and comprising 7% of the 10% total 
systematic error in the form factors. 

It is convenient to perform the chiral extrapolations at hxed pion energy, E. To do 
this, the form factors for different light quark masses are ht to the Becirevic-Kaidalov 
(BK) parametrisation [114], which satishes /+(0) = /o(0), obeys heavy quark scaling 
relations near zero recoil (g^ax) at g^ = 0 and contains a pole in /+. They are 
then interpolated to a common set of hxed values of E and extrapolated using the SyPT 
expressions [107,108]. After this, /+^o are recovered and a hnal BK ht is made to extend 
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Figure 5: Left: Fermilab-MILC-HPQCD results [109] for H —> tt and D K semileptonic 
decay form factors. The experimental points at = 0 are from BES [123]. Right: 
comparison of FNAL-MILC-HPQCD calculation [109] of /+ for D Klu with FOCUS 
data [116]. The green curve is a £t to the BK parametrisation using the lattice data. The 
green triangles show values where the fixed-U chiral extrapolation ends up. The green 
curve and points use the central values of the htted BK parameters. The blue points are 
from FOCUS. 


the results to the full kinematic range^. The output form factors are shown in hgureEland 
compared to FOCUS results [116] for Results for /+(0), compared to quenched 

calculations, are: 

Nf = 0 Nf = 2+1 
0.62(7) 0.64(3)(6) 

0.71(8) 0.73(3)(7) 

/^^^"(0)//^^^(0) 0.87(3)(9) 

The quenched resnlts are an average from [117-122], with an added 10% error incorporated 
to allow for unqnenching and lack of a continnum extrapolation. 

For decays to light vector mesons, only quenched results are available. The most 
recent results for D —> K*lu come from the SPQcdR collaboration [124]. There are four 
form factors, V and Ao,i, 2 , but Aq is not needed for the decay rate. A 2 is the least well 
determined, but since it is common to quote values at g^ = 0 , ^ 2 ( 6 ) can be determined 
from Ro(0) and Ai(0). Combining the SPQcdR results at their hner lattice spacing with 
previous results from [117,118,120,121], gives: 

U(0) = 1.17(14), Ri(0) = 0.69(7), ^ 2 ( 0 ) = 0.58(8). 

^Extrapolations have also been made using the SyPT expressions without the intermediate BK fits. 
The results are consistent but noisier [115]. 
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Figure 6: History of calculations in lattice QCD (adapted from plot shown by 
Hashimoto [138] at ICHEP04). The shaded background shows the Particle Data 
Group [22] experimental world average = (267 ± 33) MeV, with ]V[;s| as input. 

SPQcdR fit their form factors to pole/dipole ansatze (consistent with heavy quark scaling 
laws) and then integrate to find [124] (at their finer lattice spacing): 

T{D- K*^lu)/\V^^\ = 0.062(15) ps"^ 

5.1.2 Lattice Results for and fu 

Recent quenched lattice results for D-meson decay constants have concentrated on con¬ 
trolling all systematic errors save for quenching itself [130,131]. Attention has now shifted 
to simulations with dynamical quarks [126,127,139,140]. Figure IHl shows results from lat¬ 
tice calculations of published from 1998 onwards. 

Dynamical simulations with two clover flavours of light quark (although with rrisea > 
ms/2) are in progress [127]. Simulations using 2+1 flavours of improved staggered quarks 
currently allow the lightest sea quark masses to be reached. In one case [126], the heavy 
quark is treated using NRQCD and the extrapolation to the Dg meson tends to magnify 
errors. The Fermilab-MILC-HPQCD collaboration treat the charm quark directly and use 
a large set of valence and sea quark masses to help control the chiral extrapolation [125, 
140]. A ht using SyPT [107] allows discretisation errors from staggering to be removed. 
For fos, the valence quark mass is interpolated to the strange mass, and then the sea 
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quark mass is extrapolated to the up-down mass. In contrast, for the valence and sea 
masses are set equal and extrapolated together. The results are [125]: 


fo, = 249(3)(16)MeV 
Jd = 201(3)(17)MeV 


5.1.3 Sumrule Results 

Sumrules are generally less stable for D mesons than for B mesons: there are larger higher 
order operator and perturbative corrections. For the semileptonic decays of D mesons to 
TT or K mesons, Khodjamirian is updating the lightcone sumrule predictions [141] using an 
updated value for the charm quark mass, rric = 1.46 ± 0.1 GeV [22], and new Gegenbauer 
moments for the kaon and pion distribution amplitudes [142]. The preliminary results, 
compared to the summary [143] from the 2003 GKM workshop, are: 




GKM03 [143] 

2005 [141] 

f^- 

“'(0) 

0.65(11) 

0.61(11) 

/+- 

*'"(0) 

0.78(11) 

0.79(14)(8) 



0.77(4)(8) 


The second error in arises from variation in the strange quark mass, rris = 130 =F 

20 MeV. 

For the decay constants, the summary from the 2003 GKM workshop still stands [143]: 
fn = 200 ± 20 MeV, foJfD = 1-11 to 1.27 

5.2 Recent experimental results on charm decays 

Precision measurements of semileptonic charm decay rates and form factors are a princi¬ 
pal goal of the GLEO-c program at the Gornell Electron Storage Ring (GESR) [145]. We 
review herein measurements with the hrst GLEO-c data of the absolute branching frac¬ 
tions of decays to iP^e+z/e, 7r“e'''r'e ^md as well as of 77+ decays to 

TT^e'^Ue, K*^e'^Ue and p^e'^Ug, including the hrst observations and absolute branching frac¬ 
tion measurements of —*• p~e^Ve and > uje^Vf. [146,147], and compare them to 

recent measurements from other experiments and theoretical predictions. We also review 
the prospects for measuring semileptonic form factors and the GKM matrix elements Ws 
and Vcd with the full GLEO-c data set. 

The data for this analysis were collected by the GLEO-c detector at the '0(3770) 
resonance, about 40 MeV above the DD pair production threshold. A description of 
the GLEO-c detector is provided in Ref. [146] and references therein. The data sample 
consists of an integrated luminosity of 55.8 pb“^ and includes about 160,000 D^D~ and 
200,000 events. 


The technique for these measurements was hrst applied by the Mark III collabora¬ 
tion [148] at SPEAR. Candidate events are selected by reconstructing a {D~), called 
a tag, in a hadronic hnal state. The absolute branching fractions of {D^) semilep- 
tonic decays are then measured by their reconstruction in the system recoiling from the 
tag. Tagging a {D~) meson in a '0(3770) decay provides a with known 

four-momentum, allowing a semileptonic decay to be reconstructed with no kinematic 
ambiguity, even though the neutrino is undetected. 

Tagged events are selected based on two variables: AE = — Ebeam, the difference 

between the energy of the D~ tag candidate {Ed) and the beam energy (Ebeam), and the 
beam-constrained mass Mbc = \/77^eam/^"^ ~ where Pd is the measured momen¬ 

tum of the D~ candidate [149]. If multiple candidates are present in the same tag mode, 
one candidate per or D~ is chosen using AE, and the yields in each tag modes are 
obtained from hts to the Mbc distributions. The data sample comprises approximately 
60,000 (32,000) neutral (charged) tags, reconstructed in eight neutral (six charged) D 
meson hadronic hnal states [146,147], respectively. 

After a tag is identihed, the positron and a set of hadrons recoiling against the tag 
is searched for. (Muons are not used as D semileptonic decays at the '0(3770) produce 
low momentum leptons for which the CLEO-c muon identihcation is not efficient.) The 
efficiency for positron identihcation rises from about 50% at 200 MeV/c to 95% just 
above 300 MeV/c and is roughly constant thereafter, while the positron fake rate from 
charged pions and kaons is approximately 0.1% [146]. Candidates for vr® are identihd 
from photon pairs, each having an energy of at least 30 MeV, with invariant mass within 
3. Oct (a ~ 6 MeV/c^) of the known tt® mass. The Kg candidates are formed from pairs 
of oppositely-charged and vertex-constrained tracks having an invariant mass within 12 
MeV/c^ (~ 4.5cr) of the known Kg mass. K*~/p~/K *^/candidates are formed from 
{K~ and tt*^) or {Kg and vr”) / (tt” and vr*’) / {K~ and tt’*') / (vr” and combinations 
and require an invariant mass within 100 MeV/c^ and 150 MeV/c^ for the K* and p from 
the expected mean values. The reconstruction of a; —> candidates is achieved by 

combining three pions, requiring an invariant mass within 20 MeV/c^ of the known mass, 
and demanding that the charged pions do not satisfy interpretation as a K^. 

The tag and the semileptonic decay are then combined, if the event includes no tracks 
other than those of the tag and the semileptonic candidate. Semileptonic decays are 
identihed using the variable U = E-c^iss — |Pmiss|c, where E miss and pAiss are the missing 
energy and momentum of the D meson decaying semileptonically. If the decay products 
of the semileptonic decay have been correctly identihed, U is expected to be zero, since 
only a neutrino is undetected. The U distribution has a ~ 10 MeV. (The width varies 
by mode and is larger for modes with mesons.) To remove multiple candidates in each 
semileptonic mode, one combination is chosen per tag mode, based on the proximity of 
the invariant masses of the K^, K*, p, 7r°, or uj candidates to their expected masses. 

The new results for the branching ratios of various semileptonic modes are reported in 
Table El For each mode the yield is determined from a ht to its U distribution. The back¬ 
grounds are generally small and arise mostly from misreconstructed semileptonic decays 
with correctly reconstructed tags. The absolute branching fractions are then determined 
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Mode Yield B (%) B (%) (PDG) 



—>■ K e'^Ve 

1311 

.±37 

3.44 

± 

0.10 

± 

0.10 

3.58 

± 

0.18 

DO 

—>■ 7r“e’''r'e 

116.8 

± 11.2 

0.26 

± 

0.03 

± 

0.01 

0.36 

± 

0.06 

DO 

^ /C-e+i/e 

219.3 

± 15.6 

2.16 

± 

0.15 

± 

0.08 

2.15 

± 

0.35 

DO 

^ p“e+z/e 

31.1 

±6.3 

0.19 

± 

0.04 

± 

0.01 


— 


D+ 

—>• ipOg+jy^ 

545 

±24 

8.71 

± 

0.38 

± 

0.37 

6.7 

± 

0.9 

D+ 

—> TfOe+I^e 

63.0 

±8.5 

0.44 

± 

0.06 

± 

0.03 

0.31 

± 

0.15 

D+ 

^ R*0g+J,^ 

422 

±21 

5.56 

± 

0.27 

± 

0.23 

5.5 

± 

0.7 

D+ 

^ pOg+^^ 

27.4 

±5.7 

0.21 

± 

0.04 

± 

0.01 

0.25 

± 

0.10 

D+ 

—>• UJe^Ve 

7.6 

'+3.3 

'- 2.7 

0.16lO;0^±0 

.01 


— 



Table 6: Signal yields and branching fractions as recently obtained at CLEOc vs. the 
PDG results [22], The hrst uncertainty is statistical and the second systematic in the third 
column, and statistical or total in the other columns. The yield (B) for > K*~e^Ve 
is summed (averaged) over the two K*~ submodes of K~7r^ and Kg7r~. Signal yields for 

—>■ p~e^Ve and > uje^Vf. are signihcant and represent the hrst observations of 

these modes. 

using B = Ysignai/eYtag, where iVsignai is the number of fully reconstructed DD events 
obtained by htting the U distribution, is the number of events with a reconstructed 
tag, and e is the effective efficiency for detecting the semileptonic decay in an event with 
an identihed tag. 

The largest contributors to the systematic uncertainty are associated with the track¬ 
ing efficiency, the vr® and JPg reconstruction efficiencies, the extraction of D tag yields, 
the positron and hadron identihcation efficiencies, background shapes and normalizations, 
imperfect knowledge of form factors, and the simulation of hnal state radiation. The non¬ 
resonant contribution in the D —*• K*eVe reconstruction and the associated systematic 
uncertainty are accounted for as described in Ref. [146,147]. The total systematic uncer¬ 
tainty ranges from 2.8% to 7.8% according to the mode. Most systematic uncertainties 
are measured in data and will be reduced with a larger data set. 

The widths of the isospin conjugate exclusive semileptonic decay modes of the 
and are related by isospin invariance of the hadronic current. The ratio 
is expected to be unity. The world average value is 1.35 ± 0.19 [22]. Using the new 
results and the lifetimes of the and D+ [22], one obtains: \ ~ ^ 

0.05(stat) ±0.04(syst). The result is consistent with unity and with two recent less precise 
results: a measurement from BES II using the same technique [150] and an indirect 
measurement from FOGUS [151,152]. The ratios of isospin conjugate decay widths for 
other semileptonic decay modes are given in Ref. [147] and are all consistent with isospin 
invariance. 

As the data are consistent with isospin invariance, the precision of each branching 
fraction can be improved by averaging the and results for isospin conjugate pairs. 
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Decay Mode 

F (10 ^ X ps Q 

D - 

^ K e'^i'e 

8.38 ±0.20 ±0.23 

- 

->■ 7r“e^z/e 

0.68 ±0.05 ±0.02 

D - 

^ K* e+z/e 

5.32 ±0.21 ±0.20 

DO - 

p"e+z/e 

0.43 ±0.06 ±0.02 


Table 7: Isospin-averaged semileptonic decay widths with statistical and systematic nn- 
certainties. For Cabibbo-suppressed modes, the isospin average is calculated for the 
using T{D°) = 2 X T{D+). 


The isospin-averaged semileptonic decay widths, with correlations among systematic un¬ 
certainties taken into account, are given in Table [Tj 

The ratio of decay widths for D —> and D —> Ke^u provides a test of the LQCD 

charm semileptonic rate ratio prediction [153]. Using the results in Table [71 one hnds 
~ ^ 0.7(stat) ± 0.2(syst)) x 10“^, consistent with LQCD and two recent 

results [154,155]. Furthermore, the ratio ^Y(D^^e+u) predicted to be in the range 0.5 
to 1.1 (for a compilation see Ref. [151]). Using the isospin averages in Table [71 one hnds 

^rS^Ke+u) = 0-63 ± 0.03(stat) ± 0.02(syst). 

The modest in size data sample used in this analysis produced the most precise mea¬ 
surements to date of the absolute branching fractions for all modes in Table El The 
analysis already provides stringent tests of the theory. This precision is consistent with 
the expected performance of CLEO-c. CESR is running to collect a much larger "0(3770) 
data sample. 

The future CLEO-c data sample (3.0 fb“^) will allow measurements of the de¬ 
pendence (/+(g^)) and, assuming the unitarity of the CKM matrix, the absolute nor¬ 
malization (/_|_(0)) of form factors at a percent level (0(1%)) in both D —> Ke'^Ug and 
D —> Tie'^i'e [145]. For D —> K*e'^i'e and D —»■ pe'^Vei form factors A 2 {q^)., and 

U(g^) can be measured to a few percent [145]. These measurements of form factors (and 
other precision measurements in the D, B, T and 0 systems) constitute calibration and 
validation data for LQCD and other theoretical approaches. 

Theoretical predictions for form factors are required for measurements of the CKM 
matrix elements. In charm decays, the modes K~e^Ve and > 7r“e’''r'e are 

important as they are the simplest for both theory and experiment. Assuming future 
theoretical uncertainties on the form factors in these decays of 3.0% and using the antici¬ 
pated uncertainties for the absolute decay rates of > K~e^i>e and from 

CLEO-c of 1.2% and 1.5%, respectively, the following uncertainties on \Vcs\ and \Vcd\ are 
within reach: |cr(V],s)/U;s| ~ 1.6% and \(T{Vcd)lVcd\ ~ 1.7% [145]. 
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6 Conclusions 


As discussed in section El at present the determination of \Vud\ is largely dominated by 
super-allowed, O’*" —O’*" nuclear beta decays. Taking into account the new theoretical 
analysis of radiative corrections in Ref. [9], and the new global ht by Savard et al. [11], 
leads to 

\Vud\ = 0.97377 ± 0.00027 . (41) 

Hopefully, in the near future a competitive independent information on \Vud\ will be 
extracted from the neutron beta decay, once the experimental situation on the neutron 
lifetime will be clarihed. 

Concerning Vus-, at present the most reliable determination is the one obtained by 
means of Ki^, decays -see Eq. dSHI)- employing the Leutwyler-Roos value of /+(0) [49], 
which is supported by several recent lattice results [57,59-62]: 

\Vus\Ka = 0-2261 ± 0.0021 [/+(0) = 0.961 ± 0.008] . (42) 

Of comparable accuracy is the value extracted from the r(iC^2)/r(7r^2) ratio, employing 
the updated lattice result of F^/Ft^ [81]: 

|K„k„, = 0,2245;l>S Pk/F, = 1.198 ± O.OOSlSSg . (43) 

The average of these two values leads to 

decays = 0.225 ±0.001 , (44) 

which can be considered as the hnal global estimate of Vus- 

Using the values in Eqs. (pTjl and (Bill the unitarity relation in Eq. o is satished 
within la at the 10“^ level: 

A = 1 - (iKrfp ± IK/ ± IK/) = (1 ± 1) X 10-3 . (45) 
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